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Abstract 

Recent developments in the QCD description of the pion structure are reviewed. The CLEO pion- 
photon transition data analysis favors a distribution amplitude for the pion that is double-humped 
but endpoint-suppressed. After a short outline of the derivation of this amplitude from QCD sum 
rules with nonlocal condensates, we present the fully fledged analysis of the CLEO data prefaced by 
predictions for the F' 1 ^ form factor and commenting on the inherent theoretical uncertainties due 
to higher twists and NNLO perturbative corrections. We supplement our discussion by considering 
within QCD factorization theory, the electromagnetic pion form factor at NLO accuracy on one 
hand, and diffractive di-jets production on the other, comparing our predictions with the respective 
experimental data from JLab and the Fermilab E791 collaboration. In all cases, the agreement is 
impressive. 
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I. INTRODUCTION 



Advances in theoretical science are mostly based on finding useful and compact descrip- 
tions of a phenomenon of interest. Understanding the nonperturbative dynamics of quarks 
and gluons inside the pion is a question too difficult to be addressed from first principles 
within QCD. In order to gain some insight into the nonperturbative physics of the QCD 
vacuum, the concept of nonlocal condensates was introduced [1, 2, 3]. 

In this framework, a simple and compact descriptor of the complex vacuum structure is 
provided by the average virtuality (k 2 ) = A 2 of vacuum quarks, whose inverse A" 1 defines the 
correlation length of the scalar quark nonlocal condensate. The nonlocality parameter has 
been estimated using QCD sum rules [4] and lattice calculations [5]; just recently [6] it has 
been extracted from the data of the CLEO collaboration [7] on the pion-photon transition. 
One finds values in the range A 2 = (q (ig a^G^) q) /(2(qq)) = (0.35 - 0.5) GeV 2 which 
pertain to a correlation length varying, respectively, between 0.33 fm and 0.28 fm, while the 
CLEO data favors the value A 2 = 0.4 GeV 2 [6] and a correlation length of about 0.31 fm. 

Using QCD sum rules within this framework and a Gaussian ansatz with the single 
parameter A 2 for the vacuum distributions, the pion distribution amplitude (DA) was first 
derived in [1, 8]. Later, this approach was refurbished and rectified [9], providing constraints 
on the first ten moments of the pion DA and an estimate for the inverse moment using an 
independent sum rule. With recourse to the fast decrease of the moment values with its 
number N, [10, 11], the Gegenbauer coefficients of the pion DA were calculated within 
uncertainty ranges and it was found that only the first two of them a 2 and a 4 are important; 
the rest are negligible [9, 11]. This gives us a handy tool to reconstruct the pion distribution 
amplitude [9] and calculate with it pion observables, like the pion-photon transition form 
factor F ny * y *(Q 2 ) [6, 12, 13] and the electromagnetic form factor F n (Q 2 ) [14]. Confronting 
this pion DA (in terms of the first two Gegenbauer coefficients) with the constraints extracted 
from the CLEO data [7] via a best-fit analysis [15], it was found [6, 12] that it is within 
the la-error ellipse while the exclusion of the Chernyak-Zhitnitsky (CZ) [16] distribution 
amplitude at the 4cx level and of the asymptotic one at the 3a level was reinforced. 



II. NONLOCAL CONDENSATES AND THE PION DISTRIBUTION AMPLI- 
TUDE 

The pion DA is a universal process-independent characteristic of the pion, explaining how 
the longitudinal momentum P is partitioned between its two valence partons quark (x) and 
antiquark (x = 1 — x), this being a reflection of the underlying nonperturbative dynamics. 
To leading twist-2, one has 

(0 | d(z)>f>y 5 C(z, 0)u(0) | tt(P)> = if v P» f dxe ix ^<p n (x, /i 2 ) , (1) 

2-0 Jo 

where C(0, z) = V exp[—ig s f* t a A®(y) dy^] is the path-ordered phase factor (the connector 
[17]) to preserve gauge invariance, /„. is the pion decay coupling, and the normalization is 
fydxifir (x,/i 2 ) = 1. 

To get a handle on the pion DA, the aim is to relate it with the nonperturbative QCD 
vacuum in terms of nonlocal condensates. This is achieved by relating the pion and its first 
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resonance by means of a sum rule, based on the correlator of two axial currents: 

fl^) + /W*) exp = J* pZKx; s)e-'^ds + $ G (x; M 2 ) 

81M 4 ^ 1 V ' y V ; 

i=S,V,Tl, 2 ,3 

Here the index i runs over scalar (S), vector (V), and tensor (T) condensates, M 2 is the 
Borel parameter, and s° is the duality interval in the axial channel, whereas p^lo(x] s) is the 
spectral density in NLO perturbation theory [8, 9]. Above, the dependence on the crucial 
non-locality parameter A 2 enters the sum rule in the way exemplified by the numerically 
important scalar-condensate contribution 

18 ( 

$ 5 (x; M 2 ) = -j^y (x > A > x) x [x + (A - x)\n (x)] + (x -> x) + 

+6{1 > A) 9 (A>x> A) [A + (A - 2xx) ln(A)] } (3) 

with A = A 2 /(2M 2 ), A = 1 — A. One appreciates from this last expression that neglecting 
the vacuum correlation length (as in the approach of [16]), the end-point contributions 
(x — > or 1) are strongly enhanced by S(x),S'(x) . . . because for A 2 — > one obtains 
hm $ 5 (x; M 2 ) = 9 [S(x) + 5(1 — x)} . By virtue of the finiteness of A 2 , the sum rule (2) 

can supply us with constraints on the first ten moments (£ N ) n = Jq (p w (x)(2x — 1) N dx of 
the pion DA that are decreasing with increasing polynomial order to zero, i.e., (£ N ) — > 
[3/(iV + 1)(N + 3)] (see Fig. 1 (a)). In addition, we can obtain an independent sum rule to 
constraint also the inverse moment (a: -1 ),,- = J* (p^(x) x~ l dx quite accurately [9, 10] - see 
Fig. 1 (a). This is qualitatively how one derives correlated values of the first two Gegenbauer 
coefficients and in the (nonlocal) QCD sum-rules picture, given that the higher ones 
are practically zero (see Fig. 1 (b)). The final step is then to model the pion DA according 
to the expression 



^p BMS (x; pi) = <p as (x) 1 + aM) C%"(2x - 1) + a^ 2 ) Ct"(2x - 1) , (4) 



3 /2/r> i \ I /,,2\ r,3/2, 
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FIG. 1: (a) First ten (N = 10) nonzero moments, determined with nonlocal QCD sum rules [9], (£ 
and (a; -1 )^ = (1/3) (a; -1 ),,- — 1 (the superscript R meaning "reduced") of </?bms (dark blue bars) together 
with their upper and lower error-bars as light-grey bars, (b) Histogram of the first nonzero Gegenbauer 
coefficients a n of the BMS pion DA and the envelopes of the "bunch" as light- grey bars. 
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providing the "bunch" of DAs shown in Fig. 2 (Left) together with the optimum sample, 
termed BMS model [9] (02(^0) = 0-2) q 4(a*o) = — 0.14), in comparison with the CZ and 
asymptotic pion DAs. Note that all DAs mentioned are normalized at the same scale /j/q ~ 
1 GeV 2 . In Table I, we compile the main features of our SRs, contrasting them with previous 
ones. 



TABLE I: Determining the pion DA from QCD nonlocal sum rules, explaining the entries on its 
theoretical and phenomenological side, as well as its outcome. The differences between our present 
approach [9] and previous ones is pointed out. 



Ref. 


SR theor. part 


SR phenom. part 


SR estimates 


Outcome 


[1 1 


quark, quark- 
gluon and 
gluon NLCs 


7r-sidie -\- conti- 
nuum at so ~ 
0.7 GeV 2 


/e2\ /M\ lcQ\ 


r 1011 ija siiouiQ oe ciosei 
to asymptotic rather than 
to CZ DA 


[8] 


[1] + modified 
gluon NLC 


7r- and ^4i-states 
+ continuum at 
s » 2.2 GeV 2 


(£ 2 }tt,Ai, • • • , (£ 10 )7r,/li, 

(l/x)^^! with error- 
bars 


Two models of pion DAs: 
one includes the second 
Gegenbauer harmonic, i.e., 
02 7^ 0; the other is end- 
point suppressed as (xx) 2 - 3 


[9] 


[8] + corrected 
quark-gluon 
NLC Ti 


The same as in 
[8]: Borel window 
[0.5, 2] GeV 2 


(£ 2 )tt,/1i, • • • , (C 10 )7r,yli, 

(l/x) n ^A 1 with more 
conservative error-bars 


"Bunch" of self-consistent 
DAs (Fig. 2 (Left)) with 
2 Gegenbauer harmonics 
(Fig. 1(b)) 



III. LIGHT-CONE SUM-RULE PREDICTIONS FOR F 7 * p7r 

The F lp7T form factor appears as an inevitable part of the F 11 n transition form factor 
in a light-cone sum-rule (LCSR) calculation. The main advantage of this method is that 
one can calculate the form factor for sufficiently large photon virtualities to obtain the 
perturbative spectral density, and then analytically continue the result to the limit q 2 ~ 
using a dispersion relation. In this scheme, _F 7p7r expresses the "hadronic" content of the 
quasi on-shell photon 7(g 2 ) involved in the process 7*(Q 2 ) 7(<? 2 ) — > vr°. This calculational 
approach has been proposed by Khodjamirian in [18] and F lpn (Q 2 ) was computed at the 
LO level of the LCSRs. 

The form factor Q i F n ' p7V (Q 2 ) obtained in this framework depends mainly on the differ- 

d s 

ential pion characteristic —(pn(x)\ x=e , e ~ -£r, in an e-neighborhood of the origin. This 

dx Q 2 

feature is opposite to the case of the Q 2 F 11 * 7T (Q 2 ) form factor, which depends mainly on the 
inverse moment (x^ 1 }^ = ip w (x; ^)x~ x dx^ i.e., on an integral pion characteristic [9, 12]. 
From this point of view, Q F lp/K (Q 2 ) can provide complementary information on the pion 
DA and help discriminate among different pion DA models. 

Our predictions for Q^F lpn (Q 2 ), based on the BMS bunch, Fig. 2 (Left), and on a com- 
plete NLO calculation for the corresponding spectral density [6], are presented in Fig. 2 
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FIG. 2: Left:BMS "bunch" of the pion DAs contrasted with two extreme alternatives (asymptotic DA — 
dotted line and CZ model— long-dashed line) at [l 2 k 1 GcV 2 . Right: Predictions for Q i F^P 7T {Q 2 ) for the 
pion DAs shown on the left. The thickness of the two broken lines corresponds to the variation of the twist-4 
parameter in the range £>t w _ 4 = (0.15 — 0.23) GeV 2 . 

(Right) in the form of a shaded strip, with the central line denoting the BMS model. These 
calculations are sketched in the next section. Here, let us mention only the main features: 
(i) The a s -corrections appear to be rather large, of the order of 30%, and negative, (ii) The 
twist-4 contribution turns out to be very important, larger than 30% for Q 2 values below 
3 GeV 2 and also negative, (iii) An improved Breit-Wigner ansatz for the phenomenological 
spectral density p mes is used that increases the result for the form factor by about 6%. 

Comparing the different models in Fig. 2 (Left), one can understand how the slope of 
ip n (x) in the domain x ~ s p /(Q 2 + s p ) « 0.2 (at s p = 1.5 GeV 2 ,Q 2 ~ 6 GeV 2 ) translates 
into the curve for the corresponding Q i F lp ' K {Q 2 ) form factor in Fig. 2 (Right). All predic- 
tions shown are "smeared" curves, their thickness being a practical measure for the allowed 
variation of the twist-4 parameter &| w _ 4 = (0-15 — 0.23) GeV 2 . 

IV. ANALYSIS OF THE CLEO DATA 

Foremost among the many open questions in the nonperturbative regime of QCD is the 
determination of the parameters to model the shape of hadron DAs — prime examples being 
the pion Gegenbauer coefficients, on focus here, and those for the nucleon DA (for a review 
on the latter, see [19]). This task, though obviously of paramount importance is noted for 
its intransigence. However, the high-precision CLEO data [7] on the tcj transition have 
improved the situation for the pion significantly. 

Indeed, Schmedding and Yakovlev (SY) [15] have analyzed this data set using a further 
extension to the NLO (of pQCD) of the LCSR for the transition form factor F 7 * 77r (Q 2 , q 2 ~ 
0), developed in [18]. We adopted this approach in [6, 12] improving it in the following 
respects: (i) A more accurate point-to-point 2-loop ERBL [25] evolution has been employed, 
taking into account the quark thresholds, (ii) The contribution of the twist-4 term has been 
re-estimated to read <5 2 (lGeV 2 ) = (0.19 ±0.02) GeV 2 and the role of these uncertainties has 
been investigated in detail, (iii) The procedure to determine the error range of the la- and 
2a error contours has been improved and uncertainties of high-order radiative corrections 
have been involved in the analysis. For more detailed information and explicit expressions, 
the interested reader may consult [6, 12]. In the present exposition we take the opportunity 
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FIG. 3: Left: Analysis of the CLEO data on F 7r7 » 7 ((5 2 ) in the (a2, ai) plane in terms of error regions 
around the best-fit point, using for both panels the following designations: ler (solid line); 2a (dashed line); 
3tr (broken line). Various theoretical models are shown: ♦ — asymptotic DA, It — BMS model, ■ — CZ DA, 
+— best-fit point, & [20], ♦ [21], ▲ [22]— instanton models, • [23]— Ball-Zwicky (BZ) model [23], and T — 
transverse lattice result [24] . The slanted green rectangle represents the BMS "bunch" of pion DAs dictated 
by the nonlocal QCD sum rules for the value A 2 = 0.4 GeV 2 . All constraints are evaluated at fi 2 — 5.76 GeV 2 
after NLO ERBL evolution. Right: Results of the CLEO- data processing for (a; _1 )° xp /3 — 1 at /z 2 , w 1 
GeV 2 against theoretical predictions from QCD sum rules - hatched horizontal strip. The notations are the 
same as in the left panel. 

to include in our graphics (see Fig. 3) the recent results of [23] (•) and further extend our 
discussion of constraints on the Gegenbauer coefficients. 

Let us summarize the main findings of our analysis taking recourse to figures 3. (i) The 
asymptotic DA and the CZ model are excluded at 3cr and 4cx, respectively — as pointed 
out before by SY in [15]. Several other proposed model DAs, extracted, for instance, from 
instanton-based approaches [20, 21, 22], or from lattice calculations [24] are also disfavored — 
at least at the level of 2a, with the recently proposed model by Ball and Zwicky [23] lying 
exactly on the boundary of the 2a error "ellipse". The important observation here is that 
only the BMS "bunch" lies entirely inside the la error area of the CLEO data. A similar 
picture arises also for the inverse moment (x _1 )^. xp /3 — 1 (right side of Fig. 3). The light solid 
line inside the hatched band indicates the mean value of the SR estimate, (x~ 1 )^ R /3 — 1 = 
0.09, and its error bars correspondingly. The strip bounded by two almost vertical dash- 
dotted lines corresponds to the rather old Braun-Filyanov [26] constraints: <p n (l/2; fi^) = 
1.2 ± 0.3. (ii) The extracted parameters and were found to be rather sensitive to 
the strong radiative corrections and to the size of the twist-4 contribution. Nevertheless, 
even assuming a twist-4 uncertainty of the order of 30%, does not change these findings 
qualitatively. Still, both tp^y and v?cz are outside the 3a region with a slight improvement 
for the other models from 3a to 2a. (iii) The correlation length in the QCD vacuum was 
extracted directly from the CLEO data [6] and found to be A ~ 0.31 fm, i.e., A 2 < 0.4 GeV 2 . 
The prediction for the pion-photon transition form factor emerging from this analysis is 
shown in Fig. 4. The left-hand side shows the result for the twist-4 parameter S^ w _ 4 = 
0.19 GeV 2 , while the right-hand side illustrates the influence on the results of the variation 
of this parameter in the range #Tw-4 = (0-15 — 0.23) GeV 2 . As one sees, in both cases, 
the prediction for the BMS "bunch" (shaded strip) is quite robust and in good agreement 
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FIG. 4: Left: Light-cone sum-rule predictions for Q 2 F 1 * 1 - ¥V (Q 2 ) in comparison with the CELLO (dia- 
monds, [27]) and the CLEO (triangles, [7]) experimental data, evaluated with the twist-4 parameter value 
= 0-19 GeV 2 [6, 12]. The predictions correspond to selected pion DAs; notably, Lpcz (upper dashed 
line) [16], BMS- "bunch" (shaded strip) [9], two instanton-based models, viz., [20] (dotted line) and [21] 
(dash-dotted line), and <^ as (lower dashed line). Right: Similar predictions as in the left panel, but with 
the twist-4 parameter S^ w _ 4 varied in the range <5fw-4 = 0-15 — 0.23 GeV 2 . 



with both sets of data — even at relatively low Q 2 values which exceed the validity of the 
theoretical framework applied. 



V. PION'S ELECTROMAGNETIC FORM FACTOR 

We have calculated in [14] the electromagnetic pion form factor 

FAQ 2 ; /4) = F^(Q 2 ) + Fj-- WI (Q 2 ; , (5) 
where the soft part F^ D (Q 2 ) is modelled via local duality and the factorized contribution 

Fr- WI (Q 2 ; /4) = ( 2sl S +Q2 ) 2 ^ Ct (Q 2 ; /4) (6) 

with Sq~ 1oop ~ 0.6 GeV 2 has been corrected via a power-behaved pre-factor in order to 
respect the Ward identity at Q 2 = 0. In our analysis Fj act (C} 2 ; /i^) has been computed to 
NLO [28, 29], using Analytic Perturbation Theory [30, 31] and trading the running coupling 
and its powers for analytic expressions in a non-power series expansion, i.e., 

[Fl act (Q 2 ; /4)] MaxAn = «< 2) (/4) ^°(Q 2 ) + l - 4 2) (/4) ^ LO (Q 2 ; /4) , (7) 

with ai 2 ^ and (/J^) being the 2-loop analytic images of a s (Q 2 ) and (a s (Q 2 )) 2 , correspond- 
ingly (see [14] and further references cited therein), whereas T^°{Q 2 ) and J 7 ^ (Q 2 ; /x^) are 
the LO and NLO parts of the factorized form factor, respectively. The phenomenological 
upshot of this analysis is presented in Fig. 5(a), where we show F n (Q 2 ) for the BMS "bunch" 
and using the "Maximally Analytic" procedure, which improves the previously introduced 
[32] "Naive Analytic" one. This new procedure replaces the running coupling and its powers 
by their analytic versions, each with its own dispersive image [30, 31] and provides results 
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in rather good agreement with the experimental data [33, 34], given also the large errors 
of the latter. One appreciates that the form-factor predictions are only slightly larger than 
those resulting with the asymptotic DA. 
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FIG. 5: (a) Predictions for the scaled pion form factor calculated with the BMS bunch (green strip) 
encompassing nonperturbative uncertainties from nonlocal QCD sum rules [9] and renormalization scheme 
and scale ambiguities at the level of the NLO accuracy. The dashed lines inside the strip restrict the area 
of predictions accessible to the asymptotic pion DA using the "Maximally Analytic" procedure [14]. The 
experimental data are taken from [34] (diamonds) and [33] (triangles), (b): Comparison with the E791 data 
[35] on diffractive di-jet production of the BMS "bunch" (shaded strip), the asymptotic DA (solid line), and 
the CZ (dashed line) model, using the convolution approach of [36]. 



VI. CONCLUSIONS 



We have discussed pion DAs, derived from nonlocal QCD sum rules that implement the 
idea of a nonzero vacuum correlation length, and extracted for this latter quantity a value 
around ~ 0.31 fm directly from the CLEO data. The same set of experimental data offers 
the possibility to determine the Gegenbauer coefficients a 2 and within restricted error 
regions. Using a data processing, based on light-cone sum rules, we have shown that only 
the BMS-like pion DAs lie within the la region, while all other known models are excluded 
at least at 2a or more. We have given theoretical predictions for the electromagnetic form 
factor using NLO Analytic Perturbation Theory and shown that the same BMS "bunch" 
of pion DAs provides results very close to those obtained with the asymptotic DA. This 
agreement proves that whether or not the pion DA is double-humped or single-peaked is 
much less relevant compared to its endpoint behavior. In the BMS case, the root cause 
for the excellent agreement with the CLEO data is its strong endpoint suppression. These 
findings — gathered in Table II — are further backed up by the E791 data on diffractive di-jet 
production, albeit this set of data alone cannot favor one DA over the other. However, the 
fact that the middle region of x where the "bunch" has its largest uncertainties is within 
this data range, provides independent credibility for this type of pion DAs putting it into 
a wider general context. Finally, with the approved upgrade of the CEBAF accelerator at 
JLab, planned to provide precision data for the pion's electromagnetic form factor up to 
Q 2 = 6 GeV 2 , we expect that the confrontation between theory and experiment will make 
a decisive step forwards. 
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TABLE II: Pros and cons of selected pion DAs (asymptotic- like; BMS [9]; CZ [16]) in comparison 
with available experimental data. 

vr DA models asymptotic-like BMS [9] CZ [16] 
vr - 7 CLEO data [7] 3a off within la 4a off 
JLab F w data [34, 37] OK OK too large 

Fermilab E791 [35] \ 2 = 12.56 \ 2 = 10.96 \ 2 = 14.15 
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